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We consider the ground state of the homogeneous electron gas and we prove that a Hartree-Fock 
solution, motivated by previous simulations, has lower energy than the Fermi gas in the large density 
limit. This solution is a metallic phase : the density modulation corresponds to a partially occupied 
crystal (the number of sites is larger than the number of electrons) . 

I. INTRODUCTION 



The homogeneous electron gas is one of the fundamental models in condensed matter physics. Despite its simplicity 
^ - the system consists of electrons interacting with each other through a l/r potential to which a unif orm positive 
"5 background is added for charge neutrality - the phase diagram at zero temperature is nontriviaP'SEl. In general, 
^ it is given in terms of the dimensionless parameter Vg = {aodBn^ )~^, where D is the space dimension, n is the 
electronic density, as the Bohr radius and aj^ is the volume of the sphere unity. At r^ = 0, the ground state is the 
—{ Fermi gas. At large rg, the Hartree-Fock ground state of the electron gas is a Wigner crystal, that is a state where 
^ the charge density forms a triangular crystal with exactly one electron per lattice site. In the intermediate region 
$-H the Hartree-Fock approximation is not relevant and more sophisticated methods show that the ground state of the 
C/3 Hamiltonian is quite different. 

(_i At small Ts, the ground state of the electron gas within Hartree-Fock is still not known^. Although the ground 

2 state of the ideal Fermi gas (ts = 0) remains an eigenstate of the Hartree-Fock Hamiltonian of the electron gas at 
any density, already WignepSI argued that the unpolarized Fermi gas is unstable even in the limit r^ —> 0. Later, 
Overhauser showed the instability of the unpolarized Fermi gas with respect to spin-density waves within the Hartree- 
Fock approximation- . Only recently, a Hartree-Fock study of the unpolarized three-dimensional electron gas was 
performed which proposes a more complicated structure of a ground state with spin-density waves^J. 
Q In Ref .1^, we describe the results of our numerical simulations of the two-dimensional (2D) , polarized electron gas at 

'~~' small Vg, and show that the ground state is neither a Fermi gas nor a Wigner crystal: the charge density modulation 
, we find corresponds to a partially occupied crystal (the number of sites is larger than the number N of electrons) . 
^ We refer to this solution as a metallic phase. The charge density of this metallic phase represents a triangular lattice 
(^1 with reciprocal generators Qi of modulus 2kp. This modulation is mainly carried by the wave vectors close to the 
t" — Fermi surface. In order to observe the metallic phase, the number of electrons N has to exceed a threshold ranging 
r ■ from A^ ^ 10 at Ts ~ 2.6 up to A^ ^ 10^ at r^ ^ 1. As we shall see below, the threshold increases exponentially as r^ 
^^ goes to which may explain why this metallic phase has not been observed in previous simulations. 
[ — . In this study, we obtain rigorous upper bounds on the energy of the metallic phase in the two-dimensional polarized 

^^ case with a class of states mimicking the real states obtained numerically. These bounds are obtained in the limit 
00 Ts — > where the calculation is simplified by the long range behavior of the interaction potential. However, our 
numerical simulations show that this metallic phase also exists for screened potential, but strictly speaking we cannot 
!>• prove in this case the existence of such a metallic phase in the limit r^ -^ 0. This bounds are easily extended to the 
unpolarized 2D case. Then we show how this proof extends to the 3D case, whereas we have no numerical simulations 
indicating that this phase may correspond to the actual ground state of the 3D electron gas in Hartree-Fock. 



j^ Finally, let us point out that our bounds (for instance Eq. (42)) must be considered as mathematical bounds; they 

only give an understanding of the behavior of the energy of the ground state, but the actual constants involved in the 
expressions have to be evaluated by other means. 

II. NUMERICAL RESULTS 

In a previous papep^, we have computed 2D Hartree-Fock states of lower energy than the Fermi gas for values of 
Ts about 2. 

Such a state ^ is obtained as a Slater determinant Ai=i V'i where N is the number of electrons and {ipi} is a set of 
orthonormal vectors. Only the space generated by the single particle wavefunctions ipi^s is relevant, and in order to 



b(k) N=499r^=1.5 




FIG. 1: Numerical values of bk for 499 electrons in two dimensions at r^ =2. 



understand our numerical results we need to choose a canonical representation of the ipi's. 

The natural choice is to use the natural base of the Fermi gas {(j)i}i=i...N corresponding to some indexation of the 

plane waves of the Fermi sphere. 

Let M be the square matrix defined by Mij = {(pilipj). The Singular Value Decomposition (SVD) of M is M ~ UaV 

where U and V are unitary matrices and ct is a diagonal positive matrix. 

Let us define 
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(1) 

(2) 



Then {0^} and {tp'j} are also two sets of A^ orthonormal vectors satisfying {(j)i\tp'j) = Sijcr 



They are the nearest bases of Span({(/)i}) and Spand-^i}) and they define a canonical unitary operator W from 
Span({(?!)J) onto Span({V'i}) by: 



WJ2^^^^ =E.A.V^ 



(3) 



In particular, we have: 
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= T.u^3Y.^ki'^k 

j k 

= E^^^'^^fc 



(4) 
(5) 
(6) 
(7) 



Thus the natural basis of Spand-^i}) associated to the basis {4>i}i=i...N is {M^</'i}i=i...7v- 

Numerically we have chosen the ip^'s in this way and it appears that, if r^ is not too large, tpi is close to (f>i at least for 
i associated to a wave vector not too close to the Fermi surface. Thus the largest amplitude of ■0i, in the k space, is 
for the ki corresponding to (j)i. FigjT] represents the next largest amplitude of ipi, that we denote 6^., for 499 electrons 
at r^ = 2 in 2D. 

As Ts decreases, the other amplitudes of the V'i's become very small. Furthermore, the wave vector A:^ corresponding 
to 5fc. satisfies k[ — ki = Qa where the generators {Qa}a=i...6 belong to the six-fold star of a triangular lattice. 



Notice that the modulus of bk is maximal for fc on a six-fold star and close to the Fermi surface. In the direct space, 
the Qa's induce a modulation corresponding to a triangular lattice for the charge density. But while this lattice forms 
a Wigner crystal with one electron per site for r^ > 2.6, for smaller r^ it is a denser lattice corresponding to a crystal 
with an incomplete band filling. Indeed, at any Tg the electrostatic interaction favors a periodic distribution of charges 
and at small Vs only fc- vectors close to Fermi surface can be modulated. The optimal solution is to choose \\Qa\\ = '^kp 
which is larger than the reciprocal lattice vector of the Wigner case. In the following, we consider analytic solutions 
analogous to those of our numerical results and we prove that their energies are lower than the energies of the Fermi 
gas as Ts goes to 0. 

III. ENERGY OF THE FERMI GAS: POLARIZED CASE 

We consider the Hamiltonian of A^ electrons in a 2D or 3D square box of volume 17 with periodic boundary conditions. 

H^-^A+^V (8) 

2m 2 ^ ' 

where V is the 2-body Coulomb potential X^i^ii Vki ^ ''iL the electron mass is m, and e is its charge. It is convenient 
to choose Hartree as the unit of energy. Ha = fi /(ma'^), where as = ^ /{me^) is the Bohr radius. We get: 

ij = ^(-A + — y) (9) 

2 as 

Let il)n be an orthonormalized set of N vectors of _L^(fi). They define the iV-particle Slater determinant 'i! = /\„ ■0„. 
And the energy of ^ is: 

£ = (VI/|H|VI/) ^"fi-J2 {^n\A\^n) + ^ E <^" ^ ^"' 1^1^" ^ ^"') I (10) 

\ n n.n' j 

where v is defined as: 

{ipx®'^2\v\'^\®'^-il ^ dx dy lp]:{x)Tp^{y) -iji{x)tp2{y)- (11) 

In order to avoid problems due to the Coulomb singularity, we introduce the jellium model and define the potential 
acting on the plane waves 4'k as: 

^ / r, \ -D-i 

TT / Z 



(0fc(g)0fc/|w|0fe+, (g)0fe/_^) = - ( — 1 (12) 

for q ^ and otherwise, so that the total charge of the electrons is compensated by a positive background charge. 
The Fermi gas is defined by $ = A|fe|<fe 4'k where (aukp) ~ {2'K)^N/n and a^ is the volume of the unit sphere. 

\|*;|<fci=. ■" |*;|,|fc'|<fcj. ' ' 



As fl goes to oo with Cl/N fixed, the thermodynamic limit for the energy per particle is obtained by the substitution 
r 
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^k ~^ f27r)0 J ^^'- 



^^^ "^ " I I' dkk^ ^ I dkdk' ^ 



TV 2 7V(27r)^ \^7|fc|<fc^ aB2^^-i 7|fe|,|fe,|<,, \k-k'\D-^ 

" ■ "-^ ^ f dke — , \ „_^ / ^ 

|fe|<i askF^n J|fc|,|fc'|<i 



} — ^'^p^^ I / dkk^ " „ , / dkdk' ^-j—r 1 (14) 

2 Ni2n)D F IL|<, aBkF2TrD-i J,knk'\<i \k-k'\D-i ^ ' 




FIG. 2: The circle is the Fermi surface. The shaded surfaces are the regions where b{k) is nonzero. The new state tpk mixing 
(f>k and <t>k+Q^ is now resonant with t/i^/. 



From the definition of r^ = {aoaB'n} ) ^ and kp^ it follows that kpcc^fsaB = 27r. Thus, we have: 



£fg 27r2 



N 



"d ^s \"'|fe|<l 



dk k 



2 'Tsap 
47r^ 



dkdk' 



|fc|,|fc'|<i 



\k-k'\ 



J\D-l 



(15) 



which gives ior D — 2 {a^ — tt): 



£fg _ 2 

^ 7rr2 \ j|fc|<i 



dfc/c'' 



47r 



dfcdfc' 



|fc|,|fe'|<i 



l/c-fc' 



(16) 



IV. HARTREE-FOCK UPPER BOUNDS: POLARIZED 2D CASE 

We restrict ourself to the 2D polarized case and we want to estimate the energy for a class of states inspired by our 
numerical results. Let us consider a state ^ — A\k\<k ^^ where: 

V'fc = ak(j}k + bk4>k+Qt, 

with Qk in {— 2fci?(cosp7r/3,sinp7r/3)}p=o...5- For k — |/c|(cos6',sin6') we choose Qk such that |A: + Qk\ is minimal; 
that is, we choose p as the integer part of (S^/tt + 1/2) and we must assume hk is zero if k is zero or 9 = 7r/6 + n7r/3. 
Furthermore, we assume that ak and hk are real positive number and invariant thru the rotation of 2n7r/6 and the 
symmetry 9 — > —9 (i.e. the dihedral group Dq). The V'fc's are normalized, so that a^ + h\ = \ and 6/j = if 
|/c • Qk\ < 2kp{l — e) (i.e. bk is not zero only in the vicinity of {fcj?(cosp7r/3,sinp7r/3)}p=o...5), see FigJ2] 
Thus the limit energy per particle is given by: 



£2 



dk {tPkl - A\tl;k) + 



\k\<i 



in 



dkdk' —- {-tpk /\4'k'\v\')pk /\ ipk' 



|fe|,|fc'|<i 



27r 



(17) 



where, as in (16 1, the fc's have been renormalized by kp and thus \Qk\ — 2. 
We define AE by: 



£-£ 



FG 



N 



rAE 



(18) 



Then 



AE= I dk [{^k\~A\^k) 

'\k\<l 



Att 



-AEv 



(19) 



where 



AEv 



|fc|,|fe'|<i 



dkdk' i ^{tpkA V'fc'klV'fe A ipk') + t^— 



(20) 



Setting Vq = l/\q\: 



n 

2tt 



A. Potential energy contribution: AEv 



{^k A i/'fe'klV'fc A ipk') + Vk-k' = {vk-k' - Vk-k'-Q^,)bl'al 



+ {vk-k' ~ Vk+Q,-k')blal, 

+ {Vk^k- - Vk+Q,-Q^,-k')blbl, 

+ 2vQ^akbkak'bk'{6Q^+Q^, + Sq.^q^,) 

- 2vk-k'akbkak'bk'SQ^-Q^, 

— {vk+Qk-k' + Vk-k'-Qk')^kak'bkbk'6Qk+Qkf 



(21) 



Eq.([20|, may be divided into 4 parts: 

• {bk' — 0, 6fc — 0}: the contribution is zero. 

• {bk' = 0, bk ^ 0}, {bk' ^ 0, bk = 0}: both cases are equivalent. 
For {bk' = 0, fofc 7^ 0}, the integrant of Eq. 20 is: 

^i^PkAh',v^kAM + jj^^- 



{vk-k' -Vk+Qk-k')bl 



(22) 



Let So be the sector of unit disk between — 7r/6 and 7r/6 (see Fig. pi); then in this sector Qk = (—2,0) and by 
symmetry: 



dkdk' {vk-k' - Vk+Qk-k')b'i = 6 / dkdk' {vk-k' - Vk+Qk-k')bl 

bk/=0 Jk£So,bki=t) 

= 6 / dkdk' {vk-k' - Vk_k')bl 

Jk<£So,bk,=Q 

<Ce^ + 6 [ dkdk' {vk-k' - Vk_k,)bl 

JkeSo,\k'J<l-e 



(23) 
(24) 
(25) 



where fc = (2 — k^, ky). In 5*0, k — {k^, ky) where k^ is close to 1 and setting fc^: = 1 — x, we assume from now 

that bk = b{x/e). 

In Appendix A we prove that: 



L 



keSo,\k'J<l-e 



dkdk' {vk-k' - Vk+Q,-k')bl < Se'^VYe [ine^i + 0(1)] / dxb'^{x)xy/'. 



(26) 



{bk' y^O,bk^ 0}: 



By symmetry we can assume that k belongs to Sq. If k' ^ SqU S3 all the v appearing in (21 1 are uniformly 



bounded. And since the k-volume for each sector goes like e-^e, the contribution of these terms is bounded by 
Ce^. In the same way Vk-k' is bounded when k' € S3 and Vk+Q^-k' is bounded when k' £ Sq. 
Thus setting: 



/ := albl, + blal, ~ 2akak'bkbk' = {akbk' - bkak-Y 
g := albl, + b\al, + 2akak'bkbk' = (akbk' + bkak-Y 



(27) 
(28) 



one can check that: 

TT- {i'k A ipk' \v\iik A V'/c') + ^k-k' < Ce' + 6 / dkdk' {vk-k'f - Vk+Q^+k'g) (29) 

6fc,bfc,7t0 -^^ Jk,k'£So 

In Appendix B we prove that 

dkdk' {vk-k'f - Vk+Q,+k'9) < 4.e^V2~e [in e'^ + 0(1)] / dx^/i / dx' (/(ex, ex') - g{ex, ex')) (30) 

fc.fe'GSo "'0 Jx 



Thus, summing the four contribution gives: 



AE^y <Ce^ + 6e^\/2^[lne"^ + 0(l)] / dxy/x (l6b^{x)x + i j dx' {f{x,x')-g{x,x')) 



B. Kinetic energy contribution: 

The variation of the kinetic energy is given by: 







If I2 > 0, as Ts goes to 0, AE is minimal in Eq. (39) for 6 defined by 



and finally inserting ^,„in in Eq. ( 39 ) gives: 



(31) 



dk [{,Pk\ - A\,Pk) - k^] ^ & dk{{^k\ - A\i,k) ~ P) (32) 

|fe|<i JkeSo 

= 6 dx 2y„ Axb^ix/e) (33) 



< 6 X 8eV2e / dx yfx xh\x) (34) 

Jo 

C. Total energy: 

Inserting Eqs.(3l[34l in Eq. (19l, the variation of the total energy from the Fermi gas energy becomes: 



AE<6e^V2e f dx^/^ (sxb^ix) + j^ [ine-^ + 0(1)] (l6b^{x)x + 4: f dx' {f{x,x') - g{x,x')) 

= 6 X 8e2\/2^ / dxyG fxb'^{x) + ^ [ine^^ + 0(1)] fb^{x)x - a{x)b{x) f dx'a{x')bix'))jj (35) 

Let us set 

S - e^yi (36) 

/i = /" dxy/^xb^{x) (37) 

Jq 

/2 = — / dxy/^(~b^{x)x + a{x)b{x) dx'a{x')b{x')j (38) 

Then 

A£: < 6 X 8V26 [h - rshilnS-^ + 0(1))] (39) 



(Jmin^ -exp ( --^ ) (40) 



AE<-'-^^^c.J-J^]rJ, (41) 



We now have to find a solution b{x) such that I2 is positive. Choosing b{x) 
12. In the Appendix C, as Ts goes to we find a family of b leading to : 



AE < —Ts exp 



Stt 0(1) 
3r., 



60 or b{x) — 60(1 ~ x) leads to negative 



(42) 



Though such a bound is correct in the thermodynamic limit, this behavior in not so relevant for finite systems. Indeed, 
numerical systems consider about 10'^ electrons and, as we shall see later, the Fermi gas becomes then the ground 
state for r^ ^ 1. Thus the asymptotic bound (42 1 is not very helpful in the real word. 
Nevertheless, for realistic r^, on can choose a suitable function b and evaluate numerically /i and l2- For instance, 



with b = bjj as in 



of Appendix C and 77 — 0.001 we get 



A£;< -2.6 X lO^^r^exp 



18.5 



(43) 



HF UPPER BOUNDS: UNPOLARIZED 2D CASE 



We consider now the spin of the electrons, and we restrict to the case where half the electrons have a spin up and 
the others have a spin down. 
We can choose a solution as the product of 2 Slaters ^+ and ^^ for the spins up and down. 



' n n.n' 

- E (vv i^iv^n ) + — E (^" ^ V'^-'i^i^^ A ^-,) 



(44) 
(45) 



as 



E ("^n ® ^n'M^t ® ^n') + E (^" ® ^n'\^\A^ ® ^n') 



L n,n' 



(46) 



The first terms (44 1 and (45 1 can be evaluated as in (18 1 noticing that the definition of kp becomes kp = 4:nN±/il 
2TTN/fl. So the Ts in the evaluation of AE is now rsV2 and: 



£ 



NP 



cNP 



N 



AE{rW2) + m\ 



(47) 



The extra terms (46 1 coming from the direct potential provide contributions like ci^b^a^,b^^, times a positive factor. 
These contributions are regular as rg goes to 0, and does not modify the asymptotic energy. 

Nevertheless, these contributions depends on the signs of b^ in the two Slaters ^+ and ^P". Thus, since the energy 
of one Slater does not depend on the global sign of the 6^'s, the minimal state is obtained by choosing b~ = — 6+. 

For such a state, the main difference with the polarized case is that the charge density becomes flat, though the 
spin density is not. 

Finally, we can compare our solutions with previous solutions proposed by Fedders and Martin . They consider 
a more complicated modulation of the states near the center of Fermi sphere. While their results may be correct, 
they missed the dominant contribution coming from the states near the Fermi surface. Indeed, our computations give 
nonzero contribution for the states near the center, but for r^ = 1.5 bk is about 10~^ and the energy benefit, following 
their formula, must be of order 10~^" of the energy benefit of the external states. 



VI. THE METALLIC PHASE IN FINITE SYSTEMS 



For finite systems of N electrons, the minimum requirement is that the surface \k ■ Qk\ > 2kp{l ~ e), so that it 
contains at least one plane wave of the finite, discrete system. This gives the condition Ne^ > 1 and from Eqs. 
( |36|40| , this leads to: 



N > exp 



/ 3/1 
\5i2rs 



(48) 



With our approximate behavior for /2//1 in two dimensions, we find: 

N > exp ( — \ (49) 

i.e. N > 500 for r^ — 1.8. This bound is compatible with our numerical simulations where the metallic phase 
disappears at rs = 1 for N ~ 500. Furthermore, this may explain why this metallic phase has not been observed in 
previous simulations. 

VII. HARTREE-FOCK UPPER BOUNDS: POLARIZED 3D CASE 

Mutatis mutandis, the 3D case works in the same way. 
The potential term w^ is now l/||fc||^ but the dominant terms are provide now from integrals like: 



I{kx,kx')^ I dkydkzdkydk'^ —J ^^^ (50) 

||A;||,||fc'||<l 



atVyUtizOjfi.Cifi^ 11 7 /MO 



where k^ is close to k'^. 

Let us suppose that k^ > k'^, then 



Hk.,M^ Kdu)jy ,,.. , (51) 

Jo {k^-k'^y + u^ 

where u stands for {{ky — k' Y + (fc^ — fc^)^)^/^ and fi{du) stands for the distribution of u. Here we are interested in 
small values of u. One can check that /i(u, u + e) ~ 27r^(l — k^)eu provided that u is sufficiently small {{ky, k^) may 
be any point in disk of radius yj\ — fc^ and {k' k'^) may belong to a thin ring around {ky, k^) ) 

I{k,, k,,) « 7T^{1 - kl) / 2udu , (52) 

Jo {k^-k'^Y + u^ 

^-n\l~kl)\n{k^~k',f (53) 



This must be compared with the 2D case I{kx, k^') ~ — 2-\/l — k^. hi{kx — ^x)^- Mutatis mutandis, (Eq. 19 1 becomes 



AS= / dfc ((^fc|-A|^fc)-fc2)+Kr,ASy (54) 

J\k\<i 



'|fc|<i 
where n = a|/(47r^)^ and a| — 47r/3. Now, in A£^y the ^/2ex factors have to be replaced by tt^cx: 



/^Ev < Ce^ + IGzTT^e^ [ine"! + (9(1)] f dx x (b'^{x)x - a{x)b{x) f dx' a{x')h{x') ] (55) 

where z is the number of sectors of the Fermi sphere. The variation of the kinetic energy is now 

/ dk ((V-fcl - AlV'fc) - fc2) = zStte^ I dx x^b^{x) 

J\k\<l Jo 



(56) 



Thus setting 



h^ f dx x^xb^{x) (57) 

Jo 

I2 = — r— / dx X I —b {x)x + a{x)b{x) / dx'a{x')b{x') J (58) 

AE < zSire^ {h ~ ir^h [in e"^ + 0(1)] ) 



Once more, if we find a b such that I2 is positive, we have: 

AE<-—cxp(--^)rJ, (59) 

The operator A of Appendix C is now modified. But the same analysis show that the spectrum of A is still (0,4), so 
that the asymptotic bound for small r^ is: 

AE<-—exp(--^ + o{l/rs)]rs (60) 

e \ bmrrs ) 

VIII. CONCLUSION 

We have proven upper bounds for the energy of a metallic state below the Fermi gas energy. To our knowledge, 
it is the first time that rigorous upper bounds for the ground state energy of the polarized electron gas are obtained 
going beyond the simple Fermi gas states. The main interest of this proof is that it focuses on the small r^ region 
where the Hartree-Fock approximation is relevant and that the upper bounds were obtained using metallic states 
which correspond to solutions obtained by numerical simulations. 

In order to judge their relevance for the true ground state of the electron gas in the high density region, we have to 
consider correlation effects beyond the Hartree-Fock approximation. Due to the long-range, singular behavior of the 
Coulomb-potential, the perturbation expansion has to be rearranged which amounts to an effective screening of the 
electron interaction. Unfortunately a rigorous extension of our proof to the case of a screened Coulomb potential is 
not straightforward. However, numerical calculations provide evidence for the stability of the metallic states in the 
2D polarized system within Hartree-Fock, so that these states are relevant candidates for the true ground state of the 
two-dimensional electron gas. 

IX. APPENDIX A 

We have to estimate 

l-fc^<e,|*;^|<l-£ 

where k — (2 — k^, ky) and / is a positive function. 

-vr .,)= I dk'^ asinh yk^lJh. + asinh ^k±^ _ asinh ^" ~ ^^ - asinh ^" + ^^ (62) 

J f^x f^x '^x l^x '^x <^x ^ '^x f^x 

where y^ = \/l ^ ^i • And since asinh x — asinh y < Inx/y for x > y > 0: 



dk' {vu^w - vi^w) < / ' ' < 21n ^ , ^^ ,,^^ (63) 



l+e '^x k^ 



We set fca; = 1 — a; and y^ = v2a^^^2?. 



/(/) < /' dxf{x) 2y„ /' ' dk', 2 In ^^^-J^f (64) 

= 4 / dxf{x) ym / du In (65) 

Jo Je u~x 

< 4 / dxf{x) y^n I du (66) 

Jo ' Jf. u~~x 

= 8e2\/2i [In £-1 + 0(1)] / dxf{ex)xVx (67) 

Ja 



X. APPENDIX B 
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We have to estimate 

l-k^,l-k'^<e 

where k — [k^ — 2, ky) and / and g are positive functions oi \ ~ k^ and 1 — k^' ■ 
Setting kx = I — X , k'^ — 1 — x' , ky ~ y, k'y = y' , and r± = ^{x ± x')'^ + {y — v'Y, Eq 



can be rewritten: 



I{L g) = f dx y dx' I dydy' (^f- ^g\ 
= 2 I dx I dx' I dydy' i—f 



1 . 1 



where y and y' must satisfy (1 — x)'^ + y^ < 1 and (1 — x'Y + u' < 1- 
Since asinhx < ln2(a; + 1), the first term in 69 is bounded by: 



2 / da; / dx' / dydy' — / = 2 / dx / dx'/ / dy (asinh ^7 ^ + asinh ^7~^ ) 

<2 I dx I dx'f / " dy 2 asinh ^^^ 

JO Jx J -Vm X — X 

4y' 



< 4 / dx / dx'/2y„ln(2 



<4e2%/2^[ln(e-i) + 0(l)] / dx / dx f{ex,ex')y/x 

Jo Jx 



X' — X 

1 /.I 



On the other hand, using asinh x > ln2x, the last term of (69 1 is 



1 



2 dx dx' g / dydy' — — 2 / dx j dx' g / dy asinh 



> 2 / dx / dx'g / dy ln4 



X + x' 



asinh 



y'm + y 

X + x' 



/ 2 2 



And we have: 



Jx J -Vm (X + X ) 

> 4 / dx I dx'gym [hie^^ + 0(1)] 

Jo Jx 

>4e2\/2e [In £-1 + 0(1)] ( dx ( dx'g{ex,ex' 

Jo Jx 



/(/,5) <4e2^/2e[lne"l + 0(l)] / dx^/x f dx' {f{ex,ex') - g{ex,ex')) 

Jo Jx 



(68) 



(69) 



(70) 
(71) 
(72) 
(73) 

(74) 
(75) 
(76) 
(77) 

(78) 



XI. APPENDIX C 



Here we provide exact bounds on Ii and I2 given by (37 57 1 



In order to estimate I2 we introduce the hnear operator A 



Afix) 



1 



fiy)dy - 



fiy)Vydy 



^X J J. AX\J X J r\ 

defined on the Hilbert space of the functions on [0, 1] with the scalar product: 



(/l.9> = / x^/xf{x)g{x)dx 



(79) 



(80) 



11 



Then A is a bounded symmetric operator and: 



1 f {ab\Aab) 



(81) 



The unitary operator /(x) -^ g{y) — f{e ^)e ^/*^ from i^([0, l],xy/xdx) onto L'^{[0,+oo],dx) maps the operator A 
onto the operator A: 



Then 



Thus setting: 



Ag{x) = 



e/4 /-a; 



e^/4g(y)dy- 



/4 /•+00 



-y/4 



9i.y)dy 



Ae 



gk{x) 



ikx 



ikx 



-x/A 



4(1/16 + fc2) 1/2 + i2k 



1 



\l + iAk\ 
{5fe}fc>o is a full set of pseudo-eigenvectors satisfying: 



[{I + iAk)e'^'-' ~ {I - iAk)e 



■ikx~\ 



(82) 



(83) 



(84) 



Agk^ 



4(1/16 + fc2) 



9k 



Thus the spectrum of A is (0,4) and the spectral measure is purely absolutely continuous; the largest spectral value 

is 4 with a pseudo-eigenvector gi{x) = x + 4, corresponding to fi{x) — a;~^/'*(4 — lux). 

But 1 1/4 1 1 is infinite and /4 diverges at 0. The next step is to choose a family of functions fe^ such that a,, = -^1 — &,, 

is defined and (a^fe^|Aa,,5,,) /||&,)|P is close to 4. 

Thus setting fn{x) — min(/4(a;), 74(77)) for < 77 ^ 1, we have: 



Then: 



{h\Af,-,) = - 
iUAfn) 



In 7] In 1] + 0(ln 77) 

5 

3 41 2 

In rj In 7; + 0(ln 77) 

5 



l/.l 



20 

I In 7/1 



+ 0(ln-2 r/) 



Thus /,, is a good candidate for the linear part of the problem. Now, by the simple scaling: 

bjx) - -M^ 



(85) 
(86) 

(87) 



we get the nonlinear candidate satisfying 6^(a;) < l/-\/2, a^ = ^/l — bf is well defined, a^(a;) > l/-\/2 and brj satisfies 



(871 



We must now estimate the simultaneous convergence of 12/11 (81 1 and I2 as 77 decreases 



{b,,\Ab,j) - (b^aj,\Ab^ajj) = - (fe^ - fo^a^j^m - b.^a^) + 2 (fe,, - 6,,a,,|A|fe^) 

< 2(6,, - b,jan\A\br,) 

= 8 (6,, - brjUrjlbrj) + 2 (6,, - br,a,,\{A - A)brj) 

< 8 {b,, - br^a,^\br,) + 2\\b,^ - b,jaJ\\{A - 4:)b^\\ 



< 8 (6,, - br^a^lb,^) + 8||6„ - ^a„ 11^/(^1 (A - 4)6„ 



where: 



12 



{bri — brjar^\brj) — / bri{x) [1 — a,rt{x)] X\/xdx 

Jo 

f-1 



< 



' bjj{x) [1 — a,-i{x)] x^/xdxsTjp 

1 ~ Or; 



< ||6,, - 6,,a,,| 



/2-1 



and since b,j — 6^a^ > 



By direct computation: 



for 77 small enough, and thus: 



II 6,, - brfa.rfW'^ = \\bri\f - ||6r,ar,|P - 2 (6,, - 6,,a,,|6^a^) 

2 



<||6.|P 



I f^r; ^77 1 



16^11^ < II&. "' 



"" 5|ln?7| 



"■ '- -|lnr;|U/2-l5+'^', 



|6„|P 



6 l|2 



(89) 



And finally, from (87l for 6,; and (89l and h = ||&,,|p, (81 1 gives: 



OTT \ I in 77 



12 < rj^/^\ln^rj\ 

IStt 



+ 0(ln"^ r/) 
1 



1 + 



Inyyl 



(90) 
(91) 



where 



C = 20 + 8 ( /^ ^+2\/6| «92 



(92) 



Choosing ry sufficiently small, this proves that AE is strictly negative for any r^ > 0; furthermore choosing 77 to 
minimize AE (Eq. 41 1, i.e. 9rs|ln77p = 2C7r, we obtain as r^ goes to 0: 



AE < — r^ exp I — 



Stt ^ 0(1) 
3r, V^^ 



(93) 
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